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Abstract 

We derive the polarization tensor of graphene at nonzero temperature in (2+l)-dimensional 
space-time. The obtained tensor coincides with the previously known one at all Matsubara fre¬ 
quencies, but, in contrast to it, admits analytic continuation to the real frequency axis satisfying 
all physical requirements. Using the obtained representation for the polarization tensor, we devel- 
ope quantum field theoretical description for the reflectivity of graphene. The analytic asymptotic 
expressions for the reflection coefficients and reflectivities at low and high frequencies are derived 
for both independent polarizations of the electromagnetic field. The dependencies of reflectivities 
on the frequency and angle of incidence are investigated. Numerical computations using the ex¬ 
act expressions for the polarization tensor are performed and application regions for the analytic 
asymptotic results are determined. 

PACS numbers: 12.20.Ds, ll.lO.Kk, ll.lO.Wx, 73.22.Pr 


1 


I. INTRODUCTION 


During the last few years, a two-dimensional sheet of carbon atoms called graphene 


attracted much experimental and theoretical attention |l|. It is unique material whose elec¬ 
tronic excitations at small frequencies below a few eV demonstrate the linear dispersion 
relation and are described by relativistic quantum electrodynamics [l, 2|. This allows cal¬ 
culation of the reflectivity properties of graphene with respect to incident electromagnetic 
waves starting from hrst principles of quantum held theory. In fact, graphene presents the 
possibility to test at a laboratory table several effects predicted by quantum held theory, 
which were commonly considered unfeasible for observation in the past due to extreme values 
of parameters required for their realization. One could mention the interaction of graphene 
with an electrostatic potential barrier (the Klein paradox) ^ and creation of quasiparticle 


pairs in graphene either by the Schwinger mechanism in static electric held 


time-dependent held 


[0 


or m a 


d- 


Previously, the rehectivity of transverse magnetic (TM), i.e., p polarized, electromagnetic 
waves on graphene was investigated using the local model for the in-plane (longitudinal) 
graphene conductivity j?, 8|. Both the transverse magnetic and transverse electric (TE) 
rehection coefficients (i.e., for both, p and s, polarizations of the electromagnetic held) were 
expressed in terms of the polarization tensor at zero 9| and nonzero temperature. This 
approach, entirely based on quantum electrodynamics, was used to calculate the Casimir 
interaction in diherent systems involving graphene js-PZ] and turned out to be very produc¬ 


tive. Specihcally, the existence of large 


short separations, discovered in Ref. 


ihermal ehect in the Casimir force for graphene at 
18| using the longitudinal density-density correlation 


function in the random phase approximation, was conhrmed and investigated in detail. 

It should be taken into account, however, that the temperature-dependent contribution 


to the polarization tensor found in Ref. 


10| is perfectly correct only at the imaginary Mat- 


subara frequencies. The representation of Ref. 10| uses the Feynman parameter and its 
formal analytic continuation to the whole plane of complex frequency, for instance, to real 
frequencies, disagrees with that used in general thermal quantum held theory 19| . Thus, the 


results of Ref. 


10l | are well adapted to calculate the Casimir force, but are not applicable for 


the investigation of rehectivity of graphene at real frequencies. As a result, calculations of 
the rehectivity properties of graphene in Ref. 2^ using the polarization tensor was limited 
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to the case of sufficiently high frequencies, where one can omit the temperature-dependent 
part of this tensor. 

The present paper develops complete quantum held theoretical description for the rehec- 
tivity of graphene over the entire range of real frequencies from zero to inhnity. For this 
purpose we derive another representation for the polarization tensor of graphene in (2-1-1)- 
dimensional space-time which is valid not only at the imaginary Matsubara frequencies, but 
over the whole complex plane, specihcally, along the real frequency axis. Since the beginning 
of 1960s it is known which representation of the polarization tensor allows for the correct an¬ 


alytic continuation 


2l|. Here we follow methods suggested in Refs. 


HQ 


In so doing, we 


compare our representation with that obtained in Ref. 1^ and demonstrate quite a differ¬ 
ent behavior everywhere with exception of only the pure imaginary Matsubara frequencies. 
Furthermore, we consider the scattering frequency interval and hnd analytic asymptotic 
expressions for the polarization tensor and for both TM and TE reflectivities of graphene. 
The cases of low and high frequencies of incident electromagnetic waves with respect to the 
temperature are investigated. We also perform numerical computations of reflectivities and 
hnd the application regions of our asymptotic expressions. Finally, we compare quantum 
held theoretical results of this paper with those obtained previously j?, 8| using the local 
model for the conductivity of graphene. A good agreement is demonstrated for the TM 
rehectivities at both high and low frequencies. At the same time, the TE rehectivity is 
shown to draend on the angle of incidence on the contrary to a qualitative conclusion made 
previously j8|. 

The paper is organized as follows. In Sec. II we derive the polarization tensor for 
graphene applicable over the entire plane of complex frequency. Section III is devoted 
to the comparison with previously obtained expressions for the polarization tensor. The 
analytic asymptotic expressions for the polarization tensor and for both TM and TE reflec¬ 
tivities of graphene at high frequencies are obtained in Sec. IV. In Sec. V the same is done 
at low frequencies and the comparison with the results of numerical computations is pre¬ 
sented. Section VI contains our conclusions and discussion. The Appendix contains details 
of derivation of the polarization tensor in Sec. II. 

Throughout the paper we use units in which h = c = ks = ^, where ks is the Boltzmann 
constant. The fundamental constants are restored for better understanding where necessary. 
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II. THE POLARIZATION TENSOR FOR GRAPHENE AT NONZERO TEMPER¬ 
ATURE 


According to the Dirac model, the interaction of the electromagnetic field with the long 
wavelength (low frequency) electronic excitations in graphene can be described by relativistic 
quantum electrodynamics in ( 2 -|-l)-dimensions, where the ^eed of light inside the graphene 
sheet is replaced with the Fermi velocity up ~ 1/300 ll, l2|]. In this formalism, the Dirac 
equation is 

{%$ — e4 — m)ip = 0, (1) 


where m is the mass-gap parameter, is the vector-potential of the electromagnetic field, 

^ = 7^^, A = = ( 2 ) 

the Greek indices fi, u = 0, 1, 2, are the standard y-matrices and the metric tensor 77 ^ 
is given by 

= diag(l,'yF,WF). (3) 


Note that we assume graphene situated at the plane 

The polarization tensor for graphene at zero temperature was obtained in Ref. jo]. In 
Ref. [10| it was found at any nonzero temperature at all Matsubara frequencies using the 
Feynmaii parameter for a massive fermionic held and earlier, in the same representation, 
in Ref. 23(| for a massless held, but with chemical potential. Here, we obtain another 
representation for the polarization tensor of graphene which is valid over the whole plane of 
complex frequencies including the real and imaginary frequency axes. 

We start from the dehnition for the polarization tensor of graphene in the one-loop 
approximation in momentum representation lOl. l24| 


dqodq 


= STria 


tr- 


(27r)3 i(^ 


Id — m 


- 7 ' 


m 


-7 


( 4 ) 


Here, k = = (fco, k±) = [u, k±) is the 3-momentum of an external photon, ko = uj is the 

frequency, k± = {k^,k‘^), and a = ~ 1/137 in the hne structure constant. In a similar 

way, the 3-momentum of a loop electronic excitation is q = = {qo, qj_) , q j_ = (g^, g^) and 
we use the notations ^ ^ = l^k^. Note that for the purpose of calculating the 
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polarization tensor, the use of the standard 4x4 7 -niatrices instead of 2 x 2 cr-matrices and 
an additional multiple 2 included in front of the right-hand side of Eq. (jl]) take into account 
= 4 number of fermion species for graphene. 

We are interested to hnd the polarization tensor of graphene at any temperature T. In 
order to introduce temperature within the Matsubara formalism, one should replace the 
integration over go in Eq- dl]) with the summation over the imaginary fermionic Matsubara 
frequencies 

qon = mn = 27ri 0 T, (5) 

where n = 0, ±1, ±2, .... This replacement has the form 



dqo 

2tt 


00 

n=—oo 


( 6 ) 


One should also replace the frequency ko = uj in the argument of fl^'^ in Eq. (jlj) with the 
imaginary bosonic Matsubara frequencies 


koi = ui = ikii = = 27iilT, (7) 

where I is an integer number. Later it will be possible to turn the imaginary frequency back 
to real frequency. 

At both zero and nonzero temperature (as well as at nonzero chemical potential /i), the 
polarization tensor is transversal, i.e., it holds 


k^U>^%ko,k^) = 0. 


( 8 ) 


Another useful property is that in the presence of graphene all components of the polarization 
tensor can be expressed in terms of two independent quantities with corresponding form 
factors lOj. Following Ref. 10|, we use the 00-component and the trace of the polarization 
tensor 

Uooiu, k_L), Utriu, k_L) = Il^{ko,k_L), (9) 

as independent quantities in further calculations. It is only in the special case T = 0 (here 
and below we put // = 0) that the quantities ([9]) are connected by the equation 

ntr(w, k±) = - - ^ ~^^^^ noo(a;, fc^), ( 10 ) 

where k± = \k±\. As shown below, the components 
not on fc_i_ but on k±. 


of the polarization tensor depend 
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Using the Maxwell equations describing an electromagnetic field outside the surface of 
graphene and the standard matching conditions it is possible to express the reflection co¬ 
efficients on graphene for both independent polarizations via the polarization tensor. The 
result is jl^ : 


rTM(t<;, /c±) = 
rTE(t^, k±) = 

where the quantities (5 tm(te) are given by 


1 + ^-l) 

1 

1 + Qte('^) ^-l) 


k±) — noo(a;, k^), 


2ikj_ 


( 11 ) 


( 12 ) 


k±) — 


2ip{oj, k_i 

Here, we have introduced the notation 


^ . p{uj, k±) „ , , 

Htrfo;, k±) - Hoofo;, fcj 


2ik‘^, 


p{u, fcx) = 


(13) 


Substituting Eq. fll2|] in Eq. fITT]) . one obtains the reflection coefficients in terms of the 
polarization tensor [l0|, 

p(w, fc±)noo 


rTM(a;,fcx) = — 
rTE(a;, k±) = 


2ik]_ +p(a;, fc±)noo’ 

fc^Htr +P^(w, k±)Iloo 


(14) 


-2ik\p{u}, kx) + /c5.ntr + fc±)noo ’ 

where, for the sake of brevity, we put noo(tr) = noo(tr)(i^, k^_). 

Equations (IT^ and (ITTjl define the reflection coefficients (1TT|1 in the Minkowskian region, 
i.e., for real u and u > k^ such that the quantity p(a;, fcj_) defined in Eq. flT^ (it has an 
interpretation of the wave vector component perpendicular to the graphene sheet) is also 
real. The transition to the Euclidean region is done by putting u = Then Eq. flT^ turns 
into 


QTM(^e,fc±) = %|^noo(^e,fc±), 


2kl 


(15) 


<5te(*C, k±) = 




where 


2 g (^,k 

k±) = \JC + k\ = -ip{ii, kA_). 


(16) 
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Substituting Eq. flT^ in Eq. (fTT]) taken at a; = one obtains 


k±) 


rTE{i^,k±) 


q{^, fc_L)noo 

2kl + qi^,k^)noo' 

klU,,-q\^,k^)Uoo 

2g(e,fcx)fci +feint,-g2(^,fc^)noo’ 


( 17 ) 


where noo(tr) = noo(tr)(^C, k±). 

We are coming now to the calculation of polarization tensor of graphene at any tem¬ 
perature T. It is convenient to represent the independent quantities ([9]) considered at the 
imaginary Matsubara frequencies in the form 


noO(tr)(*6 , ^-l) + ^Tnoo(tr)(*6, (18) 

where A;_i_) is the 00-component (trace) of the polarization tensor at T = 0 but 

calculated at the imaginary bosonic Matsubara frequencies ife 4 ; = and A'rnoo(tr)(*6, fe±) 
is the thermal correction which explicitly depends on the temperature and goes to zero in 
the limiting case T —)■ 0. The explicit expressions for are well known 9| 


nS’Ki.U) 

U) 


ak^ 


m,k^ 


= a- 


q^{^hk_ 
q^{^i,k±) + q‘^{^i,k±) 


q‘^{^i,k±) 


m,k±), 


(19) 


where 


qi^i,k±) 

m,k±) 




m + 


2 g(C/, fe. 


Arri^ q{iu fe±) 

--arctan- 


2 m 


( 20 ) 


The expressions for A'rnoo(tr) valid at the imaginary Matsubara frequencies were obtained 
in Ref. 10|. Here, we derive the representation for A'rnoo{tr) which admits the analytic 
continuation to the real frequency axis (see Appendix A). According to Eqs. flA18p - flA2ip 
the result is 


ATnoo(tr)(*^«, fe±) 


X 



A=±l 


16a r, q± 1 

■ Jo ^^r(g^) + i 

Moo(tr){^h k±, q±) 

Ni^hk±,q±) 


( 21 ) 
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where g_L = |qj^| and the following notations are introdnced: 


N{^i,k±,q±) = [Q‘^{^i,k±,q±) - {2vFk±q±fY^^ , 

Qi^i, k±, q±) = k±) - 2iX^iT{q_i), 

r(g±) = ^ql + m^, ( 22 ) 

Mooi^i, k_i_, q±) = 4:T‘^{q±) - q‘^{iu ki_) + 4iA^/r(g_L), 

Mtriih ku g±) = -g^(6, ki_) + 4npm^ 

+ 4(l-ni) [r2(gx)+*Ae/r(gx)]. 


Note that for all ^ > 0 inclnding the Matsnbara freqnencies with I > 1 the snmmation in 
A = ±1 in Eq. (|2T]) can be eqnivalently represented as 


1 il4oo(tr) 

2 ^ N 

A=±l 


Re 


^OO(tr) 

N 


(23) 


We emphasize that in order to obtain the correct resnlt in Eq. (l2l]) at = 0, one should first 
make the summation in A = ±1 at .^ > 0 and then consider the limit ^ 0. Note that if one 

puts ^0 = 0 hrst, the subsequent summation in A = ±1 would lead to a complex result for 
some relationships between k± and q± in contradiction with the fact that the polarization 
tensor along the imaginary frequency axis must be real. The proposed continuation of the 
expression (EB to zero Matsnbara frequency is equivalent to the use of Eq. fl23|) at all 
.^ > 0. It is also conhrmed in Sec. Ill by the comparison with another representation for the 
polarization tensor of graphene. 

In the massless case, m = 0, the above representation for the polarization tensor of 
graphene can be simplihed. Using Eqs. flT^ and (1201) . for the zero temperature parts one 
obtains 


fzU k<) — 

T-r{0),.^ , ^ q^{.ihki_) + q'^{iuk^) 

Ill'{i^i,k±) = na - -. 

k±) 

For the thermal correction in the massless case T{qj_) = qj_ and from Eq. 


(24) 


we have 


N{^,,kug±) = «&,n)4r(6,n,«±), 


(25) 


where 


W{ii, i±, ?±) = [?^(^i, ki) - iiXiiqL - iqi]■ 


( 26 ) 







In a similar way, from Eq. fl2^ in the massless case one arrives at 


Moote, q±) = k^, gx), (27) 

MU^i, fcx, q±) = -(1 - vl)W\^i, fcx, q±) - fcx). 


Substituting Eqs. fl25l) and fl271) in Eq. fl2Tl) . we find the thermal correction to the polarization 
tensor in the massless case 


^Tnoo(*65 k±) — 


16a 


dq± 


Jo 


ei±/T + 1 


X 


. _ ^ k±, q±) 

. «teXi) 


ATntr(f^i, k±) = 


16a 


dq± 


-'f -'0 


ggx/T _j_ I 


X 




A=±l L 


(1 - 


q{^h k±) 


+ 1 ) 1 ; 


qj^i, k±) 


(28) 


Finally, we consider the Minkowskian region by turning the imaginary frequency to the 
real frequency axis in accordance to a; = Here we have several frequency subintervals. 
For the lowest frequencies, u < all the above formulas hold literally and q remains 

real after putting The interval where surface plasmons may exist goes next. It is 

characterized by 


upfex < OJ < fcx, (29) 

q{u, ki_) = + k\, p{uj, fcx) = \J- vlk\, 

where q and p are real. The scattering interval 


kj^ <oj 


(30) 


is the most important one because it corresponds to real photons on a mass-shell. Here, 
both the quantity p(a;, /cx) dehned in Eq. flT^ and p(a;, fcx) dehned in Eq. fl29|l are real. 

In the plasmonic and in the scattering frequency intervals, the polarization tensor has a 
similar form. At zero temperature it is given by 


nS’(wXi) 


ak\ 


<I)(a;, fcj 


= a- 


p 2 (x;, ki) 

p^(a;,/cx) +P^(cu, fc±) 

p\jj,kx) 


<I)(a;, /cj 


(31) 


9 

















where the function <h along the real frequency axis is dehned as 


/c_l) 


1 L - arctanh%^’ 


m 






to < ^/vpk‘j_ + Am?, 

^ _ (32) 

u > \/v^k\ + Am?. 


The thermal correction to the polarization tensor in the plasmonic and scattering frequency 
intervals takes the form of Eq. fl2T]) . 


A j'nQQ(tr) /c_l) 


16« r , q± 1 

vl Jo ^^V{q^)en<i.)IT + l 


X 


^ 1 ^ 00 (tr)(l^, <?±) 

2 ^^^ N{uj,k^,qi_) 


(33) 


Here, the denominator N takes different forms depending on the value of the quantity 
[compare with Eq. fl22|) ] 


Q{u, /c_L, g_L) = -p^{uj, k±) - 2XuT{q±). (34) 

If \Q{uj, q±)\ > 2vFk±qj_ we have 


N{uj, k_i_, q±) = signQ(a;, k±, q±) 

X [Q‘^{uj, k±, q±) - {2vFk±q±Y] , (35) 

whereas if \Q{u;,k±,q±)\ < 2vFk±q± 

N{u, k±, q±) = -i [-Q‘^{u, k±, q±) + {2vFk±q±f] (36) 

holds. Two other quantities entering Eq. fl55]l are given by 

Moo{uj, k±, q±) = 4r^(g_L) k±) + 4An;r(gx), 

k±, q±) = k±) + 4npm^ (37) 

+ 4(l-n2) [r2(g^) + Aa;r(g^)]. 


Note that in the plasmonic frequency interval the denominator N has no zeroes, but in the 
scattering interval with A = — 1 it has two zeros at 


{q±)i,2 — T 


VFk± 

2 




Am? 

p 2 (a;,fc^)' 


These zeros result in the integrable singularities. 


( 38 ) 
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Now we consider the massless case in the scattering interval where Eqs. fl33|) and (l35|) -(l37 |) 
for the polarization tensor can be signihcantly simplihed. At zero temperatnre the resnlt 
follows from Eqs. fl3ip and fl32p in the limiting case m —)■ 0 20l |. 




mak^ 


kA_) ’ 

p2(a;, fc_L) +p'^{u:,kA_) 


-moL- 


p{uj, kx 


(39) 


For the thermal correction in the massless case from Eqs. fl35l) and fl36l) one obtains 


N{uj,k^,qx) = p{u,k_L)W{uj,kx_,qi_), 


(40) 


where 

W^{u, k±, q±) = p‘^{u, k±) + 4Aa;g_L + 4g^. (41) 

The explicit expression for W{oj,k±,q±) follows from Eqs. (I3^ and fl36p . It is different for 
different A = ±1 and depends on the valne of q±. Thns, for A = 1 and for any q± it holds 

IE(a;, k±, q±) = - [p^{u, kx_) + 4a;g_L + 4g^] , (42) 


whereas for A = — 1 we have 

- [p‘^{uj, k±) - 4:uq± + 4:ql ]^^‘^, 

-i [-p^iu, kx_) + 4a;g_L - , 

k±) - 4uq± + , 


W{u,k±,q±) = < 


2q± < uj — VFk±, 


oj — vpk± < 2g_i_ < a; + vpk±, 


2qx_ > a; + vpk^. 


(43) 


In so doing all expressions under the square roots in Eqs. fl4^ and fl43P are positive. 

Using Eq. gB, the quantities Mqo and Mtr dehned in Eq. fl37H can be written in the 
massless case as 


Moo{uj, k±, q±) = IU^(a;, k_i_, q±), (44) 

Mtr{u, fc_L, q±) = (1 - fc_L, q±) + vIp‘^{u}, kx_). 

Substituting Eqs. fl40P and fl44p in Eq. ([33]), we arrive at the thermal correction to the 
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polarization tensor of graphene with m = 0 in the scattering interval, 


16a 1 

ATlloo{u,k^) = ^ dq^ 


JO 


eq±/T + 1 


X 


^ 1 ^ W (g;, fc_L, q±) 


( 45 ) 


^ ^ , 16a /■“ , 1 

ATntr(a;,/cx) = —2“ / “?-L 


7)2 / 

% Jo 


X 




A=±l L 


(1 - ^f) 


+ 1 

hh(w, fc L, g L) 


p(w, fcj 


+ np 


P(^) fe±) 

W{oj,k^,qA_) 


Note also that in the scattering interval (1401) under consideration here, the denominator of 
A-rfltr has zeros at 


{q. 


uj ± npfcx 


-Ln ,2 - 


(46) 


These singularities are integrable. 


III. COMPARISON BETWEEN DIFFERENT REPRESENTATIONS FOR THE 
POLARIZATION TENSOR OF GRAPHENE 


In previous section, we have obtained the polarization tensor of graphene at the imaginary 
Matsubara frequencies and found its analytic continuation to the whole plane of complex 
frequencies including the real frequency axis. As discussed in Sec. I, another representation 
for the polarization tensor is contained in Ref. [^. Here we compare both representations 
and demonstrate that the tensor of Ref. [ij can be used solely at the Matsubara frequencies, 
but not along the real frequency axis [i.e., Eqs. (13) and (24) in Ref. Q]. For the sake of 
brevity, here we restrict ourselves to the case of gapless graphene (m = 0), which is used in 
our calculations of reflectivities in the following sections. 

First of all we note that the polarization tensor of graphene at zero temperature calculated 
at the imaginary Matsubara frequencies is given by Fq. (|24)) in both representations. Because 
of this, the subjects of our comparison are only the thermal corrections A'rnoo(tr)- Taking 
into account the structure of the reflection coefficients dm), below we compare different 
representations not for Apnoo(*6 , A;_l) and ATnti.(i^i, fc±) but for Apnoo(*Ci, A;_l) and the 
following quantity: 


Apn(i^,, ki_) = kA.) - q {iu k±)ATUoo{i^i, k±). 


(47) 


12 















In our representation A 7 ’noo(*^n ^J.) is given by the first line of Eq. 
expression for ATn(z^i, k±) follows from Eq. 


AxUii^i, k±) = 


16a 


-'f Jo 


dq± 


ggx/T _|_ I 


-ei 


and the explicit 


(48) 


..2 7^2 q{^hk 


+q te,fc±)Re- 


m,k±) 


We also take into account Eq. fl^. Our representation should be compared with the fol¬ 


lowing alternative results of Ref. 


10l | (see also Ref. which we notate by an upper index 


(a): 


ATn;,“^(*en k^) = ^ r dx\T\n[l + 2 cos(e;x/T)e-®^("’«*’"^i 


-^F Jo 


, -20x(x,^„fcx)l _ _ 2r) _ M^iJ^/T) _ 

^ cosh 0(r(a;, /cj_) + cos(.^/a;/T) 

— x) cos(^/x/T)-I- 1 
q{^hk±) cosh0r(x,^«,fc_L) + cos(^/x/T) J ’ 


(49) 




8 a 


-^F -'0 


dx \ T^f In [1 + 2 cos(Cza:/T)e-®^(^’^*’^^i 
i^i/4 sin(^ix/T) 


where 


^g-2exfe5„x)]_[2.j2(5,,A^J_^2]2(i_2i)—— 

2 cosh 07 ’(a;, 4;, fc_i_)-|-cos(4;a;/T) 

rrt-; 4 cos(4ia;/T)-4 I 

+ Vx{l- x)q^{^i, k±) - c i. - (c IT\ I ’ 

cosh Qt[x, 4/, /c_l) + cos[iix T) j 


OT(x,0,k±) = ^x(l - x). 


(50) 


As shown in Refs. 


10 


12 j, for graphene the major contribution to the thermal correction 


comes from the zero Matsubara frequency. We start from this case and by putting 1 = 0, 
4 o = 0 in Eqs. (1281) and (I48|) obtain 

dq± ^ r^Fk ±/2 /^,2 u 2 _ A ^2 


« ^ /„ , N 16a 

lkTnoo( 0 , k±) = —^ 


A'rn( 0 , k±) = 16a 


'0 

f*vpkj_/2 


gijx/r _|_ 1 vpk± Jq 
dq_i ( vpk^ 


1 - 4,3 


ggx/r _|_ I 


(51) 


ki 


+ 1 I ^vlkl - 4g2 


+ —Jvlkl-Aql 

v-p ^ 


By introducing the new variable ^ = 2g_i_/(nFfc_i_), these results can be rewritten in the form 


ATnoo(0, fcx) = fTln2 - ^ C dz-'^ 


gA 2/2 _|_ ]_ 


A(rn( 0 , fcj_) = — 8 anFfci 


dz 


+ 1 yrW72 ’ 


(52) 


13 



































where A = v-pkj^/T. 

The alternative expressions of Ref. [l^ are obtained by putting Z = 0, .^o = 0 in Eq. 

8 aT 


ATn|,j>(o,n) = 


dx In 


T -'0 


Arn*'“^(0, fcj_) = —/ dx\/x{l — x) 




( 53 ) 


JO cosh[74y'a:(l — x)] + 1 

By making transformations in the hrst line of Eq. fl53l) and integrating by parts, we have 

16aT 


ATnS( 0 ,fcx) = 


dx In 


T JO 


X _|_ g-Ay'ajll-x) 


IbaT 




T\n2 + VFk± f ^ dx 


a;(l - 2x) 

+vpk± / dx 


0 2y^x{l^^ gAi/xp-a:) 

1 


h /2 2\/x{\ — x) gA-^x(l-a:) ^ 

In both these integrals we introduce the new variable 

l-2x 


(54) 


y 


= 2^/x{l - x), dy = 


a/x (1 - x)’ 


Xl,2 = 


i± 


y" 


(55) 


where the signs plus and minus correspond to the second and hrst integrals, respectively. 
This leads to the result coinciding with the hrst line of Eq. fl52|) . 

In a similar way, the second line of Eq. fl5^ can be represented in the form 


A'rn^“^(0 , A;_l) = —IGavpk^ [ dx- ^ 


'0 e 


A^x{l-x) _|_ 


= —Ibanp/c^ 


^Jx{l - x) \/x{l 

dx ^ r + / dx- ^ ^ 


X 


(56) 


/o gAy/x(l-x) y ^/2 gA\/a;(l-x) 

After introducing the new variable (l55ll . where the sign plus corresponds to the second 


integral and the sign minus to the hrst, Eq. fl56|l coincides with the second line in Eq. 
Thus, it holds 

ATn„(0, k^) = Arn‘;>(0,4x) Arn(0, fci) = ATn(“>(0, k^). 


(57) 


In fact our representation for the polarization tensor in Eqs. fl28|) and fl48|) coincides with 


that of Ref. 


10| [see Eq. fl4^ ] not only at zeroth but at all Matsubara frequencies. To 
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illustrate this, in the second column of Table I we present the values of ArllgQ^, normalized 


by the coefficient C = 16aT/u|, which are computed using the first formula of Eq. flT9|) at 
fcx = 10.^1 and first eleven nonzero Matsubara frequencies indicated in the first column. In 
the third column of this table, our results for Arlloo computed using the first formula of 
Eq. (j28il are presented. As can be seen in Table I, both sets of results coincide up to a high 
accuracy. The fourth column contains the relative thermal correction AT^IIoo/ngQ^ which 
decreases quickly with increasing 1. Similar computations were performed at = O.S^i, 
100.^1 and 500.^i. It was found that in all cases the values of Arnoo (*65 computed using 
different representations coincide. Note that with respect to the Casimir effect 25|, l2^ the 
value k±_ = lO^i is the magnitude of the wave vector giving the major contribution to the 
force between two parallel graphene sheets spaced at the separation distance 100 nm. 

In Table II similar computational results for the quantities k±) and 

are presented normalized by the coefficient D = Sa^f/vp. As can be seen 
in this table, the thermal correction At-II computed using the representation of Ref. 

[i.e., the second formula of Eq. (14^ . column 2] and the approach of this paper, i.e., Eq. (l48ll 
(column 3) coincide at all Matsubara frequencies (column 1). The same values of k±, as in 
Table I, were used in computations. The relative thermal correction is given in column 4. 
Thus, one can conclude that 


kj_) = ATTl‘m\i(h U) A 3 .n(!{,, kj_) = A 3 .n(“>(ij,, kj_). 


(58) 


It can be easily seen, however, that the polarization tensors of graphene obtained in 
Ref. |l^ and in this paper behave quite differently at all frequencies other than the Mat¬ 
subara frequencies. To see this, in Fig. 1 we plot 10^ x At^oo/C (a) and AtIIqq^/C' (b) as 
functions of at k±_ = lO^i. In so doing, the quantities Arlloo and Arllgg^ were computed 
using the first formulas of Eqs. fl28l) and fl49l) . respectively, where the Matsubara frequencies 
are replaced with the continuous frequency ^ [see Eqs. (32) and (33) of Ref. |27| where 


the polarization tensor of Ref. 


10 | is written along continuous pure imaginary frequencies]. 


As is seen in Fig. l(a,b), our thermal correction decreases monotonously with increasing 
whereas the thermal correction of Ref. 10| oscillates between the values it takes at the Mat¬ 
subara frequencies. This is a nonphysical analytic continuation from the set of Matsubara 
frequencies to the imaginary frequency axis and, as a consequence, to the whole plane of 
complex frequencies including the real frequency axis (see Sec. IV for the physical results). 
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In Fig. 2(a,b) we present the respective results for the quantity ATn(i^, fcj_) at fcj, = lO^i. 
Computations were performed in the framework of our formalism by Eq. fHHjl [Fig. 2(a)] and 
using the formalism of Ref. by the second formula of Eq. 0491) [Fig. 2(b)]. The quantities 
10^ X A'^n/D and Aj'Il^^^ /D are again plotted as functions of continuous frequency along the 
imaginary frequency axis. In our formalism, dependence of the thermal correction A-pIl on 
the frequency is again monotonous, whereas the quantity oscillates taking the same 

values as AtH only at the Matsubara frequencies. The respective continuation of ATn*^“) 
to the real frequency axis does not possess necessary physical properties (see Sec. IV). This 
is an example of the well known general situation in the complex analysis that a sequence 
/(/) with no accumulation point can give rise to different analytic functions F{z) such that 
F{1) = f{l). For example, the function F(z + 1) provides an analytic continuation of n\ to 
the entire complex plane, but this continuation is not unique. In particular, the analytic 
function, e. g., of the form F( 2 : + 1) + (7 sin(7rz), where C = const, also provides the desired 
continuation. 


IV. REFLECTIVITY OF GRAPHENE AT HIGH FREQUENCIES 


Here, we hud analytic asymptotic expressions for the polarization tensor of graphene in 
the scattering frequency interval fl5UD and respective graphene reflectivities at high frequen¬ 
cies satisfying the condition ca S> T. If we restore the fundamental constants, this condition 
takes the form 


(jj t/Jx = 


keT 


h 


(59) 


where ut is so called thermal frequency. For example, at room temperature T = 300 K the 
thermal frequency ut ~ 3.9 x 10^^ rad/s, i.e., Eq. fl59l) is well applicable in the optical range 
and at all higher frequencies. 

Calculations are done starting from Eq. fl45l) for ATnoo(w,/c±) along the real frequency 
axis and the following expression for A 7 ’n(ci;, fc^), entering the reflection coefficient tte in 
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Eq. (fT^ . which is obtained from Eq. fH5|) : 


Arn(a;, /c_l) = k\/\T^tr{oo, k^) + p^(a;, A;_L)A'rnoo(a;, fcj 


16q! 


dg L 


ggx/r + 1 


UJ 


p‘^{u:,kx) ^ fE(a;,fc_L,9±) v^k^ ^ 


A=±l 


A=±l 


hh(w,/c_L, g_L) 


( 60 ) 


At zero temperature the quantity IIqo^ is given by the hrst line in Eq. (13^ and the quantity 
n(o) is given by 


U^^\u,k^) = klU^^,\u,k^) +p\u,k^)uS}{u,k 


( 0 )/ 


= —i7rak‘]_p{u, k±). 


(61) 


We note that at /c_l = 0 the quantity Arnoo(iA;,0) = 0. This is seen from Eqs. (021), 
(03|1 and (051) if to take into account that at A = 1 iy(a;, 0, gj,) = —oj — 2q± whereas at 
A = —1 kE(a;,0,g_L) = 2q± —u holds for any q±. Then, by using Eq. (1601) . one concludes 
that A'rn(a;,0) = 0 as well. As can be seen from Eqs. fl39D and (l6T|) . the same property 
nQg^(a;, 0) = n(°)(a;, 0) = 0 holds at zero temperature as well. 

Now we obtain the asymptotic expression for At^IIoo under the condition fl59D . The main 
contribution to the integral in the hrst formula of Eq. (051) is given by q± ~ T. Taking into 
account that u: ^ T and u: > k±, one can conclude that the second and third intervals of 
q_[_ in Eq. 03i) do not contribute to the result essentially. Keeping in mind that for g_L ~ T 

OJ — 2q± ^ OJ ^ VFk± (62) 


holds, for A = 1 from Eq. 02|l we obtain 


Similarly, from the hrst line of Eq. 03ll we hnd 

2 7 2 

W_i = lV(u,k±,qJ ^ -(u- 2q^) + 

2(oj - 2q_L) 

Taking into account that 

p(cj, k±) OJ \ 2oj^ / ’ 


(63) 


(64) 


(65) 
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( 66 ) 


from Eqs. 


Substituting Eq. 
part of A-rlloo: 


one arrives at 

1 fEi + fE_i 


-1 + 


2vlklql 


2 p{u, k±) 

in the first formula of Eq. fH5|) . one obtains the hnal result for the real 


ReATlloo{u;,k±) ^32a^ dq± 

^ Jo 




eq±/T + 1 


L 2 rp 3 

= 48aC(3)^, 




(67) 


where ((’(z) is the Riemann zeta function. 

In the above, we have neglected by the imaginary part of ImArlloo which originates 
from the second line of Eq. fH3|) . It can be shown, however, that under the condition fl5^ 
ImAThloo is exponentially small. Really, from the second line of Eq. fl43|l and the hrst 
formula of Eq. fl45p we have 


ImATnoo(tu, k±) 


8a ri^+v^k^)/2 (2g^_^)2]V2 

vjpiu, k±) + 1 


( 68 ) 


Introducing the new variable t = 2q±/{vFk±), this can be represented in the form 

, , ^ , 4afc2 /-/hFfcx)+i 

ImATnoo(tu, k±) = 


p{u, kj_) e^Ffcxi/(2r) + 1 


X Wl- U- 


u 


-V 


VFk± J 


(69) 


Due to Eq. ([62]), this is an integration over the narrow interval centered at u/{vFk±). At 
the ends of this interval the square root in the integrand is equal to zero and in the middle 
to unity. Thus, the trapezium estimation of the integral results in 


ImATnoo(a;, k±) ^ 


Aak^ 1 

- v^kl e-/(2D + 1 


(70) 


conhrming an exponential smallness of the imaginary part of the polarization tensor at high 
frequencies. 

We are coming to calculation of the reflection coefficients and reflectivities of graphene 
in the scattering region (l30ll . For photons on a mass-shell the condition k^ + k^ = uP' is 
satished, so that fc_L can be expressed as 


k^ = oj sin 0*, 


(71) 
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where 6i is the angle of incidence. From Eqs. fl3^ and fl671) we obtain the following expression 
for the 00-component of the polarization tensor at high frequencies 


noo(i^, kj 


iTrak'^t k'^tT^ 

+ 48aC(3)- 


(72) 


p{u, k±) “ ' ' 

Substituting this in Eq. flT^ for the TM reflection coefficient and using Eq. fTflll and the 
inequalities vp 1, T/u ^ 1, one obtains 

Tia cos 6i(l + ^ sin^ 0^ — 48ia( (3) ^ cos 




2 + 'K a cos 6i [l + ^ sin^ 6i ) — ASiaC (3) ^ cos 6i 


Tia cos 9i 


,^^F.2n 48C(3)/T 

1 H-sm Oi — i - ' 

2 TT 


UJ 


(73) 


2 -|- 7ra cos 6i 

The reflectivity of the TM polarized light from graphene under the condition (15^ follows 
from Eq. fl73|l 


|r tm(i^, 9i)Y 


COS^ Qn 


1 I 2-2^1 /48(C(3)^ 
1 -|- Up sm -|- I -j 


(74) 


(2 -|- 7ra cos^j)^ 

Note that the corrections to unity in the square brackets of Eq. fl74)) are very minor. This 
is seen from the fact that Up = (up/c)^ ~ 10“® and T/o; = ksTj^hx)) ~ 10“^ at room 
temperature and optical frequencies. Because of this, in the high frequency region the 
relative thermal correction to the TM reflectivity of graphene at room temperature is of 
order 10“^^ (at zero temperature the result (1741) was obtained in Ref. (20j). 

In Fig. 3, we plot the TM reflectivity of graphene computed by Eq. fl74)) as a function of 
Qi (the lower line). At the normal incidence (6*j = 0) we have 


|?"tm(i^, 0)1' 

in agreement with the result of Ref. 


TT^CK^ 


(2 -|- ttq;)^ 


1.28 X 10 


-4 


(76) 


8| obtained numerically using the local model for the 
conductivity of graphene. As is seen in Fig. 3, at high frequencies the TM reflectivity of 
graphene is rather small and decreases with increasing angle of incidence. 

To consider the TE reflectivity of graphene at high frequencies, we need the asymptotic 
expression for the quantity Apll defined in Eq. fl60|) . Substituting Eqs. (I63|) and (16T)) in 
Eq. fl60|) . one obtains 


ReA(rn(a;, k±) ~ —32akj 


qI 


^±/T + 1 _ g2 


(76) 


-32a^ 


=-48aC(3)qA. 


e<i±/T + 1 


■a;" 
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Similar to the case of ImArlloo, it is easily seen that ImATll is exponentially small. 

As a result, the asymptotic expression for If at high frequencies is obtained from Eqs. fl6T|) 
and (17^ 

T3 

n(a;, kj_) ~ /cj_) — A8a({3)k\—. (77) 

We substitute this equation in the second line of Eq. flT^ and use Eq. (ITT!) and the smallness 
of our parameters up and T/ca: 


r te(i^, di) 


Tra 1 - ^ 


Tiaoj — ^sin^^jj — 48iaC(3)^ 

2u cos 9i + TTauj (^1 — ^ sin^ 9^ — 48iaC (3) ^ 

f r _ . 48C(3) /tV' 


) 




(78) 


The reflectivity of the TE polarized light from graphene at high frequencies is given by 

2 /m\ 6 ' 


\rTEiuJ,9i)\‘^ 




(2 cos 9i + na) 


1 2-2 
1 — Up sm 




(79) 


Similar to Eq. (|7411 . the corrections to unity in the square brackets of Eq. (17^ are very small. 
At T = 0 the result fl79l) was derived in Ref. 20]. 

The TE reflectivity of graphene computed by Eq. (ITOj) as a function of 9i is shown in 
Fig. 3 (the upper line). At the normal incidence {9^ = 0) we have 


|?"te(i^, 0)1^ — |?"tm(i^, 0)1^, (80) 

where |rTM(i ^5 0)p is given in Eq. fl75|) . This is in agreement with Ref. js]. However, with in¬ 
creasing angle of incidence, |rTE(i^, 9i)\'^ monotonously increases, as opposite to |rTM(i^, 
which decreases with the increase of 9i. Thu^ at this point our quantitative results dis¬ 
agree with the qualitative conclusion of Ref. [8| that in the framework of the Dirac model 
\rTE{9i)\^ = |r’TM(0)P holds for any angle of incidence. 


V. REFLECTIVITY OF GRAPHENE AT LOW AND INTERMEDIATE FRE¬ 
QUENCIES 

In this section, we derive analytic asymptotic expressions for the polarization tensor and 
reflectivities of graphene in the scattering frequency interval (130|) at low frequencies satisfying 
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the condition u: <^T. We also perform numerical computations using the exact expressions 
for the polarization tensor and hnd explicitly the application regions for the asymptotic 
expressions at low and high frequencies. 

We start from the imaginary part of the thermal correction At^IIoo in the hrst formula of 
Eq. fl45p , which is determined by the second line of Eq. (145]) 


^ ^ , 8a riu^+vyk ^)/2 

ImArllooiu, k±) = — — —— / „ . 

v^p{uj, k±_) + 1 


8a 


vIp{u, k±) 

+ / 


ruj/2 




{uj - 2q^Y (81) 


dq± 


■\Jvlk\ - (cj - 2qi_y 


{uj-v^k^:)/2 + 1 

{u)+VYk^_)/2 1 . - 


lui/2 


It is convenient to introduce new variables x = u — 2qj_ and x = 2q± — ca in the hrst and 
second integrals on the right-hand side of Eq. (IHTj) . respectively. This results in 


ImArnoo(a;,/ c_l) 


4a 

k±) 



dx 


.2U2 


2 + e 2 t -\- e~ 


x^ 


iJ — X 

e 2 t 


+ 1 


e 2 t 


+ 1 


(82) 


Now we take into account that x uj because in the scattering region under consideration 
VFk± k± < u. In the lowest order one can put u ±x ^ u and using the condition cj -C T 
obtain 


ImATnoo(w, fcj 




—vra- 


(83) 


p{uj, k±) e2T -1- 1 

Comparing this result with Eq. (139|) . we conclude that at low frequencies the thermal correc¬ 
tion to IIoo gives the same contribution to the imaginary part of IIoo as the zero-temperature 
part n^. 

The real part of the thermal correction At^IIoo is determined by Eq. fl42]) and by the hrst 
and third lines of Eq. f|43l) . Substituting these in the hrst formula of Eq. (1451) . one hnds 


16a 

ReATnoo(w, fc_L) = 


dq± 


gijx/T _|_ 1 2p{u), k±) 


dq_L 


\/{u + 2qi y - v^k 


2p{u, kx 


A(a; - 2qx_f - v^k 
dq± — 


2p{uj,k^)^ 

I {uj+vpk±)/2 


dq± 


ei±/T + 1 

2qx_)‘^ - vjkl 

ggx/T _j_ I 


(84) 


In this equation, we make the following identical transformations. The hrst integral is 
calculated explicitly. In the second and third integrals we introduce the new variables 


21 






























y = oj At 2gj_ and y = uj — 2gj_, respectively. In the fourth integral we introduce the new 
variable y = 2q± — u. Then Eq. fl8T|) takes the form 


16a 

ReATlloo{u,k±) = 


-'F L 


Tln2- 


4:p{u, k±) + 1 


^ _. A,2 _ „.2 

/(2r) 1 V ^ 


1 


dy 


Ap{uj, k±) e^oj-y)l{2T) ^ X 

1 r dy 


Ap{u, fcx) Apfcx + 1 


^1/2 - vlk\ 


(85) 


In the last integral on the right-hand side of Eq. fl8^ let us separate the integration region 
oo) into two parts: [vFk±,u;] and [a;, cx)). Then we combine the integral over the hrst 
part with the second integral on the right-hand side of Eq. fl85l) and the integral over the 
second part — with the hrst. As a result, we have 


16a r 1 /- 

ReArn„((.,,fcJ = ^ T\n2 + —— dy^y^-v^^kl 

^F L ^±) Jv^kj^ ^ 

e-‘^/(2r) ^y/{2T) 

_j_ g—a;/(2T) ^yjip.T') ^Lojip.T') 

4p(w, fc±) L V ^ l^e2^/(2^) -K 6-/(2^) eS'/(2F) + e—/(2r) 

This is an exact result which can be used both for numerical computations and for obtaining 
the asymptotic expressions at low frequency. 

It is easily seen that at a; T the hrst integral on the right-hand side of Eq. fl86|l 
is negligibly small. To see this we expand the quantity in round brackets in porews of 
small parameters u/T and y/T and hnd the leading order term equal to —y/{AT). After 
integration with respect to y the magnitude of the hrst integral in Eq. fl86l) together with 
its coefficient appears to be v‘^k\/{^AS)T). 

The second integral on the right-hand side of Eq. fl86|) can be calculated taking into 
account that w -C T and ^ ^ < y- Then we can expand in the small parameters 

VFk±/y and u/T. Keeping in mind also that the main contribution to the integral is given 
by 1 / ~ T, we can replace the lower integration limit with zero. Then from Eq. fl86p we 
obtain 





ReApnoo(w, k±) 


16a 


Tln2- 


Aoj 



roo Qy/{2T) 'j 

[ + i]2 j 



2a;2 J 


(87) 


22 




























and after integration by parts arrive at 


ReAT’noo(i^, ~ —8aTln2 






( 88 ) 


Thns, the contribntion of the second integral in Eq. fl86|) is really ranch larger than that of 
the Erst. 

As a resnlt, from Eqs. fl8^ and fl88l) one obtains the following expression for the 

00-component of the polarization tensor at low freqnency 


noo(tt’, — 




A’2 

8aTln2^. 




(89) 


ATp{uj, k±) 

Now we snbstitnte Eq. flS^ in Eq. ffn)) . use Eq. flTT]) and replace p{u,k±) with u. In the 
case of low frequencies this replacement is appropriate because the terms of order are 
much smaller than the thermal correction. Then we have 


6*j) 


Tia-^ cos 9i + Aia In 2 ^ cos 6i 
1 -|- TTO)^ cos 6i + Aia In 2 ^ cos di 


(90) 


From Eq. (I90|) the reflectivity of the TM polarized light from graphene at low frequencies is 
given by 

'jr‘^ — cos^ O' “t" 16 q)^ ln^2 —— cos^ O' 

|rTM(a^.9i)P« d , “L • . (91) 

(l -|- ttq;^ cos -|- lOa^ ln^2 ^ cos^ Oi 

Note that at low frequency both real and imaginary parts of the IIoo component flH^ con¬ 
tribute to the reflection coefficient fl9UD and reffectivity flTID . The thermal correction con¬ 
tributes essentially to both quantities. At the normal incidence we have 


|rTM((^, 0)1" 


+ 16a^ ln"^2 


2m2o 


,2o 7E 
.,2 


(92) 


(l -|- TTQ!^) -|-16Q!2 1n^2 

In Fig. 4, we plot the reflectivity of graphene (19^ at = 0 as a function of frequency 
measured in K at three different temperatures T = 10 K, 100 K, and 300 K (the lower, middle, 
and upper lines, respectively). As is seen in Fig. 4, at a; = 0 the reffectivity is equal to unity 
and decreases monotonously with increasing frequency. The comparison with the results of 
numerical computations (see below) shows that even at the lowest considered temperature 
T = lOK the asymptotic expression (19^ is rather precise. Our analytic results of Fig 4 
obtained using the polarization tensor are in rather good agreement with the computational 
results of Ref. (see their Fig. 2, left) found using the local model for the conductivity 
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of graphene. Note, however, that with decreasing temperature the local model supposedly 
overestimates the TM reflectivity of graphene. 

We are coming now to the calculation of the TE reflectivity of graphene at low frequency. 
Besides floo, the reflection coefficient tte in Eq. flT^ is expressed via the quantity fl dehned 
in Eqs. (1^ and fET]) . To calculate Aj’!! at a; <C T, it is convenient to use its dehnition in 
the hrst line of Eq. fl60|) and already obtained results fl83|) and (l88|) for AtHoo- The quantity 
Arfltr is given in the second formula of Eq. fHSjl . Both its imaginary and real parts at low 
frequency are calculated using the same procedure, as was applied above in this section to 
calculate ImArfloo and ReA-rfloo- The results are 


ImA'rntr(c(;, fcx) ~ tto ' 
ReA7’ntr(a.;, fcx) ~ 8Q!Tln2 




p{uj,k±) e^ + l’ 

k±) +p‘^{u, k±) + 2vpk'j 


(93) 




Then from the second line of Eq. 


fltr ~ 8aTln2 


and Eq. fl93|) one hnds 
p^(a;, k±) +p^(a;, k±) + 2vlk‘j_ 




.u p^{u,k±) +p^{u,k±) 

+ * 777 ;^“-——-— -. 

4T p{uj, k±) 


Finally, from Eqs. (18^ and (l9T|l we arrive at 


(94) 


n(a;, fcx) ^ 8aT ln2 k‘]_ — i—ITak‘]_p{u, k±), 


(95) 


where we have taken into account that Vpk"^ <C 

Substituting Eq. fl95|l in Eq. flTT)) for the TE reflection coefficient, using Eq. flTB and 
taking into account that p{uj, k±) ^ a; in the scattering frequency interval, one obtains 


rTE{(^,0i) = 


TTCK^ + 4mln2 - 

_ OI _L(J_ 

cos R + TTCK^ + 4*0; ln2 - 

01 (jJ 


(96) 


The reflectivity of the TE polarized light from graphene at low frequencies is given by 


|rTE(i^, 6'j)p — 


+ 16q;^ ln^2 ^ 


(cos^i + TTo;^) + 16 q; 2 ln^2 ^ 

From the comparison of Eq. (177)1 with Eq. (19^ it is seen that at the normal incidence 


(97) 


|?"te(i^, 0)1^ — |’"tm(i^, 0)1^, 


(98) 
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as it should be, where the right-hand side of this equation is given by Eq. ([92]). Thus, 
Fig. 4 demonstrating dependence of the reflectivity of graphene on frequency at the normal 
incidence is equally applicable to both TM and TE polarized light. 

Now we illustrate dependence of the reflectivities of graphene at low frequency on the 
angle of incidence. Calculations were performed using Eqs. fl9ip and fl97p for the TM and 
TE polarized light, respectively. In Fig. 5 the calculation results are presented by the upper 
pair of solid lines computed at u/T = 0.01 and lower pair corresponding to u/T = 0.03. As 
is seen from Fig. 5, at low frequencies the reflectivities of TM polarized light decrease with 
increasing 0*, whereas the reflectivities of TE polarized light are the monotonously increasing 
functions of 0*. In this sense the dependence of reflectivities of graphene at low frequencies 
on the angle of incidence is similar to that at high frequencies. The single difference is that 
the magnitudes of reflectivity at low frequencies are several orders of magnitude higher than 
at high frequencies. 

In the end of this section we present the results of numerical computations for the re¬ 
flectivity of graphene over a wide range of frequencies and compare them with the analytic 
asymptotic results obtained in this and previous sections. As an example, we choose the in¬ 
cidence angle 9i = 30° and perform numerical computations of the TM and TE reflectivities 
as functions of cu/T in the wide range from 10“^ to 10^. Let us start from the TM polarized 
light. In this case computations are done using Eqs. dn, (USD, ini) and fl45p . The com¬ 
putational results are shown by the solid line in Fig. 6(a) plotted in the double logarithmic 
scale. In the same hgure, the long-dashed line shows the analytic asymptotic results at high 
frequencies. They are calculated using Eq. (1741) at 6^ = 30°. The analytic asymptotic results 
at low frequencies are calculated by Eq. (1^ and shown by the short-dashed line. As is seen 
in Fig. 6(a), the analytic asymptotic results are in a very good agreement with the results of 
numerical computations over a wide frequency region. Thus, at low frequencies the relative 
deviation of the short-dashed line less than 5% occurs at all u/T < 0.3 [log(ci;/T) < —0.52]. 
At high frequencies less than 5% relative deviation holds at all u/T > 9 [log(a;/T) > 0.95]. 

Finally we consider the TE polarized light incident on graphene at the same angle 6i = 
30°. We perform numerical computations using Eqs. ffT4|) . (1431) . fl60l) and ([6T|). In Fig. 6(b), 
the computational results are presented by the solid line plotted in the double logarithmic 
scale. Similar to Fig. 6(a), the analytic results of Eq. d79|l at high frequencies are shown by 
the long-dashed line and the analytic results of Eq. (1971) at low frequencies are shown by 
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the short-dashed line. From Fig. 6(b) it is seen that at low freqnencies less than 5% relative 
deviation occurs at u/T < 0.25 [log(a;/T) < —0.6], and at high frequencies less than 5% 
relative deviation holds at all co/T > 9 [log(a;/T) > 0.95]. This solves the question when one 
can use simple asymptotic formulas and when more cumbersome numerical computations 
are necessary. 


VI. CONCLUSIONS AND DISCUSSION 


In the foregoing, we have developed quantum held theoretical description for the rehectiv- 
ity of graphene described in the framework of the Dirac model. This description exploits the 
polarization tensor of graphene at nonzero temperature in (2-|-l)-dimensional space-time. 
The previously known representation of this tensor obtained in thermal quantum held 
theory in Matsubara formulation is applicable only at the imaginary Matsubara frequencies 
and does not admit physically reasonable analytic continuation to the real frequency axis. 
Here we derive another representation for the polarization tensor of graphene which pos¬ 
sesses the required analytic properties over the entire plane of complex frequencies. At zero 


9j . At nonzero 
10] at all Matsub 


temperature our polarization tensor coincides with that obtained in Re 
temperature our tensor takes the same values as the one obtained in Ref. 
ara frequencies. This justihes all numerous applications of the polarization tensor of Ref. [K 
to calculation oUhe Castottt .ntecaction iu ayatema tucludiug gtaphene and gtaphene-coapd 
substrates [ICHlTl]. However, as shown in this paper, the polarization tensor of Ref. [1 


is an oscillating function along the imaginary frequency axis, and this precludes the use of 
its analytic continuation to real frequency axis for theoretical description of the rehectivity 
of graphene. On the contrary, in the representation obtained here the polarization tensor 
is a monotonously decreasing function of the imaginary frequency and admits the analytic 
continuation to real frequency axis satisfying all physical requirements. 

Using the derived representation for the polarization tensor at nonzero temperature, we 
have obtained analytic expressions for the reflection coefficients and reffectivities of graphene 
for two independent polarizations of the electromagnetic held in the asymptotic regions of 
high and low frequencies, as compared to temperature. For the transverse magnetic {p 


oolarized) Ih 


it our results are in a good agreement with those previously obtained at high 


8(] and low [7| frequencies using the local model for the conductivity of graphene. However, 
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some qualitative conclusions about the reflectivity of TE (s polarized) light made previously 
are not supported by our exact calculations. We have investigated dependences of the TM 


and TE reflectivities of graphene as functions of frequency and incidence angle. In the latter 


case it was shown that the TM reflectivity is a decreasing function of the incidence angle, 
whereas the TE reflectivity increases with increase of the angle of incidence. 

Furthermore, we have performed numerical computations of the TM and TE reflectivities 
of graphene over the wide frequency range. By comparing the analytic and computational re¬ 
sults the application regions for the asymptotic expansions at low and high frequencies were 
determined. It was shown that simple analytic expressions for the reflectivity of graphene 
obtained here produce very accurate results in the wide regions of frequencies and temper¬ 
atures. 

By and large the developed formalism solves the problem of the reflectivity of pristine 
(gapless) graphene and graphene characterized by some nonzero mass-gap parameter. In 
future it would be interesting to generalize this formalism to the case of nonzero chemical 
potential and consider applications to doped graphene sheets. 

Appendix A: Polarization tensor in (2-|-l)-dimensions at nonzero temperature 

In this Appendix, we display the details for the calculation of the polarization tensor of 
graphene at nonzero temperature. By making the replacement (E]) in Eq. (jlj), we represent 
the polarization tensor in the form 



OO 



n=—oo 


where the spinor propagator is given by 



— m — m? — iO 


(A 2 ) 


and the notations introduced in Eqs. ([5]) and ([ 6 ]) are used here and below. 
By calculating the trace over the 7 -matrices, we get 



(A3) 
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where 


q) = VX' [/(/ - fc"') + (/ - ~k^')q^' 

-g„(g“ +mV'l , 

R{k, q) = {qf,q^ - + iO) (g^ - k^,){q^ -k^)-m‘^ + iO 


(A4) 


In there equations, = {iki 4 ,VFk±) and is the metric tensor in (2+l)-dimensions with 
= —g^^ = 1. Note also that we have made the substitution q_i_ —)■ q^/up in the 


integration in Eq. flAip . 

For the independent quantities (jO]) characterizing the polarization tensor at any temper¬ 
ature it follows 


noO(tr)(^) 


327ra f dc[± Zoo(tr)(A:, g) 

vl (27r)2 R{k,q) 


(AS) 


Here, by using Eq. flA4l) . we obtain 


Zm{k, q) = -g4n(g4n - kii) + q±(q± - Upfe^) + 
Ztr{k, q) = -(1 - 2np)g4„(g4„ - k^) + q_L(q_L - 
+ (1 + 2n|)m^, 

R{q. k) = [gL + r^(?±)] \{qln - kh) + T‘^{k±, q±)l , 


(A 6 ) 


where T{q±) is dehned in Eq. (j2^ and 

f (fc±, q±) = [(q± - Vwkxf + m^] 


(A7) 


We mention that Hoo in 


coincides up to an overall factor with Eq. (A3) in Ref. 281. 


Now, following the idea behind the Abel-Plana formula (see, e.g.. Ref. 25(]), we represent 
the Matsubara sum over n in Eq. (IASI) as an integral in the complex g 4 plane, 

327ra f dq^ f dq± 1 Zoo(tr)(fc,g) 


n 


OO(tr) 


{k) = 


2ti 


71 U 72 


(27r)2 1 i?(/c, g) 


(AS) 


where the integration path 71 U 72 encircles the poles of the real axis at the fermionic 
Matsubara frequencies indicated by dots (see Fig. [7]). In the same figure, the zeros of the 
quantity R(g, k) in the denominator of flASp at g 4 = iiE and g 4 = ±iT + k^ are shown as 


crosses. 






















It is convenient to represent Eq. (lASp as the sum of two integrals. The hrst of them is 
along the path 71 and the second is along the path 72 . In order to have a decreasing function 
in the upper half plane under the integral along the path 71 , we use the identity 


1 

gig4/T _|_ I 


1 - 


1 

Q-m/T X 


(A9) 


We substitute Eq. flA9p in the hrst integral of Eq. (lASp along the path 71 . In so doing the 
unity results in the polarization tensor at zero temperature 


(0) . . ^ 327ra f f dq_L ^oo(tr)(fc,g) 

Vp J 2ti J (27r)2 R(k, q) 


(AlO) 


_ _32to ^ f dq_L ^oo(tr)(fc,g) 

27r y (27r)2 R{k,q) 

The substitution of the second term of Eq. flA9D in the hrst integral of Eq. flASp along the 
path 7 i together with the second integral of Eq. flASp along the path 72 result in the thermal 
correction to the polarization tensor 


. T-r ^2'Ka f dq^ f dqi_ 1 Zoo(tr)(A:, g) 

^ Vp J 271 J ( 27 r )2 e“*'?4/T _|_ i R(^k, q) 

71 

327ra C dq^ C dq± 1 Zoo(tr)(/c, g) 

Vp J 271 J (27r)2 0 * 94 /'^ + 1 R(^k, q) 

72 


(All) 


It is well known that the polarization tensor has no ultraviolet divergencies in (2+1)- 
dimensions. In fact, its divergence degree is +1, which would correspond to a linear di¬ 
vergence. The latter, however, is not present due to the gauge invariance. Otherwise, the 
linear divergence would reside in the zero-temperature part IIqqI^j.^. The thermal corrections 
A'rnoo(tr) do not have any ultraviolet divergencies, which is evident due to the Boltzmann 
factors. 

Calculation of the polarization tensor at zero temperature (lAlOp was performed in 
Ref. fITop with the result flT^ . Now we calculate the thermal correction to it (lAlip . For 
this purpose, we close the integration pathes 71 and 72 by the half circles of inhnitely large 
radii in the upper and lower half planes, respectively. These circles do not contribute to the 
result which is given by the corresponding residua in the poles of both integrands shown in 
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Fig. [7] as crosses 


A TT 327ra 

Arnoo(tr)(fc) = —— 


ciq_L 

( 2 vr )2 


E 

A=±l 


er(qx)/r + i 


^oo(tr)(fc,*Ar(g L),q L) 


2Ar(g^)[(zAr(g^) - + V^{k^,ci^)] 

1 Eoo(tr)(fc, ^4;+ *Af(fc_L,q_L),q±) 


(A12) 


g[_iAfc4,+r(fc^,qx)/T + 1 2Ar(fcx, q±) [{iXT{k^, q^) + k^i)^ + F 2(g^)] 

where the contributions with A = 1 come from the upper half plane and those with A = — 1 
from the lower half plane. 

Eouation flA12D can be simolihed takine into account that according to Eo. o ^ 4 / are 
the Matsubara frequencies and, thus, exp(—iAfc 4 ;/T) = 1. In fact doing just this way is 
necessary for a correct analytic continuation to non-Matsubara frequencies since otherwise 
the replacement k4i —)■ ioj would not result in a function flA12l) decreasing for oj —)■ ±cxd 
as is required for the analytic continuation. In application to temperature Green functions 
this was hrst mentioned in Ref. (2l|] and extensively elaborated in Ref. 22|. A recent, 

[see Eq. (A5) there]. Further 


mathematically equivalent, treatment can be found in Ref. 


we use the symmetry of the integrand under the substitution q_i_ —)■ ~ q± and making 

the replacement A —)■ —A in the second term of Eq. flA12ll we arrive at 


ATlIoo(tr)(fc) — 


IbTTtt f dq± 1 


(27r)2 r(gx) + 1 


(A13) 


X 


E 

A=±l 


^oo(tr)(fcGAF(gx),q±) + Eoo(to)(fc, ^4/ - ^AF(g_L),- q±) 
[k^i - iXT{q±)]‘^ + F^(fc±,q±) 


As said above, in this representation, the analytic continuation of Apnoo(tr) from the Mat¬ 
subara frequencies k^i to any frequency k^ = ca, including the real frequency axis, can be 
done directly. 

Using the polar coordinates (gj_ = |qj_| = ql + q^-, </?) in the qj_-plane, it is possible to 
perform one of the two integrations in Eq. flAlSp . Taking into account that 


2 qj_fcj_ = 2gj_fcj_ cos(p. 


(A14) 


for the quantity in the denominator of flAlSp we obtain 


[kii - iAF(g_L)]^ + F^(fcx,q±) = <5(^40 - Iv^q^k^cosip 


(A15) 
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with the notation 

Q{hi, k±, q±) = kli + vlk^ - 2i\kiiT{qi_). 

In the numerator of Eq. IA13P we get 

Zoo{tr){k,iXT{q±),q±) + Zoo(^tr){k, k^i - iXT{q±),VFk± - q^) 
= Q{kii,k±,q±) - 2vFq±k±cos qp + Moo^tr){kAi, k±, q±) 


(A16) 


(A17) 


with 


Mooikii, k±, q±) = -(kh + vlkl) +AiXk4iT{q±) +4r^(g_L), 

Mtrik^i, k±, q±) = -{kh + vlk]_) + 4npm^ (A18) 

+ 4(1 - vl)iXk4ir{q_F) + 4(1 - n|)r^(g_L)- 

Rewriting the integral flAlSp in the polar coordinates and using Eqs. flAlSD and flA17D . 


we obtain 


, ^ ,,, 167ra 

^Tnoo(tr)(^) 2 / 


dq± 


(1± 


X 


p 2 tt 

A=±l “'O 


(i(p 


T{q±) er( 9 x)/r + 1 
Moo(^tv){k4i, k±, q±) 


(A19) 


0 271 Q{k^i,k±,q±) - 2vFq±k_i_ cos If 


The remaining angular integral can be carried out by using the formula 

dip 1 1 


Jq 277 Q - a cos p - a- 

with a = 2vFq±k±. Here, the branch of the square root must be dehned such that the sign 
of the imaginary part of the root be the same as that of Q. This prescription 

holds for Q as given by Eq. flAlbp with real k^i. For a complex k^ one needs to use the 
corresponding analytic continuation. Substituting Eq. flA20p in Eq. flAlOp and using the 
notation 

N{kii, /c_L, q±) = [Q‘^{kiu q±) - {2,VFq±k±fY^‘^, (A21) 

one arrives to the hnal result (IMil for the polarization tensor of graphene dehned at the imag¬ 
inary Matsubara frequencies. The analytic continuation of this tensor to the real frequency 
axis is described in Sec. II. 
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TABLE 1: The normalized values of the thermal correction to the 00-component of polarization 
tensor of graphene computed at the Matsubara frequencies with different numbers (column 1) using 
the representation of Ref. {l^ (column 2) and in this paper (column 3) at = lO^i. The relative 


thermal correction is given in column 4. 
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Arnoo/n™ (%) 
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0.24 

4 

3.6872960104 x lO"'^ 

3.6872960109 x lO"'^ 

0.108 

5 

1.547382154 x 10“’^ 

1.547382155 x 10“’^ 

0.056 
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7.56717034 x lO"® 

7.56717035 x 10"® 
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4.12036077 x 10“® 

4.12036078 x 10"® 
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1.52276096 x 10“® 

0.010 

10 

1.00201160 X 10“® 
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6.85884190 x 10“® 
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TABLE II: The normalized values of the thermal correction (BZI) to the combination of the com¬ 
ponents of polarization tensor of graphene computed at the Matsubara frequencies with different 
numbers (column 1) using the representation of Ref. [l^ (column 2) and in this paper (column 3) 


at k± = 10 ^ 1 . 

The relative thermal correction 

is given in column 4. 


i 

ATn(“)/D 

AtH/D 

Ayn/ni®) (%) 

1 

1.12028612025 x 10“"^ 

1.12028612024 x lO'^ 

4.08 

2 

4.02785688588 x 10“® 

4.02785688597 x 10"® 

0.73 

3 

2.0002626999 x 10“® 

2.0002626997 x 10“® 

0.24 

4 

1.1797934387 x 10“® 

1.1797934384 x 10“® 

0.108 

5 

7.736289303 x 10"® 

7.736289297 x 10“® 

0.056 

6 

5.44805008 x 10“® 

5.44805006 x 10“® 

0.033 

7 

4.03778027 x 10"® 

4.03778024 x 10"® 

0.021 

8 

3.10943615 X 10“® 

3.10943614 X 10“® 

0.014 

9 

2.466807364 x 10"® 

2.466807362 x 10“® 

0.010 

10 

2.00397985 x 10"® 

2.00397988 x 10"® 

0.0073 

11 

1.65980994 x 10"® 

1.65980992 x 10"® 

0.0055 
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FIG. 1: (Color online) The normalized values of the thermal correction to the 00-component 
polarization tensor of graphene computed as functions of imaginary frequency at k±_ = lO^i 


this paper and (b) using the representation of Ref 


. Q. 


of the 
(a) in 
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FIG. 2: (Color online) The normalized 
the components of polarization tensor 
at k_\_ = 10^1 (a) in this paper and (b) 


values of the thermal correction (1471) to the combination of 
of graphene computed as functions of imaginary frequency 


using the representation of Ref, 


. Q. 
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FIG. 3: Reflectivities of the TM and TE polarized light of high frequencies on graphene are shown 
as functions of the incidence angle by the lower and upper lines, respectively. 
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FIG. 4: (Color online) Reflectivities of graphene at the normal incidence at different temperatures 
T = lOK, 100 K, and 300 K are shown as functions of frequency measured in K by the solid lines 
from bottom to top, respectively. 
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FIG. 5: (Color online) Reflectivities of the TM and TE polarized light of low freqnencies (w/T = 
0.01 and oj/T = 0.03 for the upper and lower pairs of lines, respectively) on graphene are shown 
as functions of the incidence angle. 
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FIG. 6: (Color online) Reflectivities of the (a) TM and (b) TE polarized light on graphene are 
shown as functions of uj/T at the incidence angle Oi = 30°. The solid, short-dashed and long- 
dashed lines are computed by the exact formulas, and by the asymptotic expressions for low and 
high frequencies, respectively. 
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FIG. 7: The complex ( 74 -plane containing the integration pathes 71 and 72 . The dots indicate the 
poles at the fermionic Matsubara frequencies and the crosses indicate four additional poles (see 
text for further discussion). 
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